Memory Effects and Scaling Properties of Traffic Flows 
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Traffic flows are studied in terms of their noise of sound, which is an easily accessible experimental 
quantity. The sound noise data is studied making use of scaling properties of wavelet transforms 
and Hurst exponents are extracted. The scaling behavior is used to characterize the traffic flows in 
terms of scaling properties of the memory function in Mori-Lee stochastic differential equations. The 
results obtained provides for a new theoretical as well as experimental framework to characterize 
the large-time behavior of traffic flows. The present paper outlines the procedure by making use 
of one-lane computer simulations as well as sound-data measurements from a real two-lane traffic 
flow. We flnd the presence of conventional diffusion as well as 1// noise in real traffic flows at large 
time scales. 
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I. INTRODUCTION 



Various aspects on traffic flows has been studied in 
great detail over the years (see e.g. Refs.[ll2 and ref- 
erences cited therein) . Due to the properties of vehicles 
and the variability of drivers behavior, it appears quite 
natural that stochastic methods is the proper theoreti- 
cal framework in modelling of traffic flows. Numerical 
. simulations based on the Nagel-Schreckenberg scenario 

I— ' I d, 3 have provided us with various insights in the large 
scale behavior of traffic flows and various actual analyt- 
ical modelling has been discussed in great detail the lit- 
erature (for some recent accounts see e.g. Refs. 5]). 

It has become an increasing demand in populated ar- 
eas to monitor and control the sound noise level of traffic 
flows. It is an easy experimental task to collect sound 
data of a traffic flow, but it is our understanding that 
such data has not until now been used to study the prop- 
erties of the traffic flow by itself. In the present paper we 
introduce what we believe is a new and novel experimen- 
tal method to study some characteristic features of traffic 
flows which can be combined with the theory of stochas- 
tic differential equations of the Mori-Lee type 0, Q with 
memory effects included. In the course of our discus- 
sion below we will see that this actually can be a fruitful 
method. By making use of a now well established scal- 
ing analysis in terms of wavelets, one can infer self-affine 
properties of traffic flows. Generic properties of Mori- 
Lee 0, stochastic differential equations are then in fact 
compatible with the scaling properties obtained from the 
data analysis. We use artificial traffic data as obtained 
from one-lane computer simulations with car stopping in- 
cluded as well as data from a real field measurement of 
an opposite two-lane traffic flow with a common cross- 



road in order to test our procedure. The real traffic data 
exhibits short-time scaling corresponding to normal Ein- 
stein diffusion, a finite time 1// noise (see e.g. ReL 8] and 
references therein) behavior as well as large-time fluctu- 
ations of the sound-level corresponding to a flat power 
spectrum. 

The paper is organized as follows. In Section^]we re- 
call some basic properties of stochastic differential equa- 
tions of the Mori-Lee form in terms of conventional 
Laplace transform techniques. The large-time properties 
of fluctuations in the noise level are expressed in terms 
of the short distance properties of the Laplace transform 
of the memory function. In Section IlIII we analyze time- 
series of sound measurements of traffic fiows. The scaling 
properties of the traffic data is based on a now estab- 
lished method of wavelet scaling properties of statistical 
data and the extraction of the so called Hurst scaling 
exponent 9, 10, llj. In Section llVI we summarize our 
results and suggest further extensions of our analysis of 
traffic fiows. 

II. MEMORY EFFECTS AND THE MORI-LEE 
EQUATION 



We consider the dynamics of the traffic system in terms 
of the sound noise they produce. This noise will be mod- 
elled by a second-order stochastic differential equation 
for the time-dependent sound intensity /(<) treated as 
a stochastic variable with a space-homogeneous memory 
function of the form [12j | 



dv{t) 
dt 



+ / dt'r{t - t')v{t') = F{t) 



(1) 
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Here v{t) — dl{t)/dt/l{t), and F{t) is a stochastic back- 
ground. Eq.(^ can be derived from the wave-equation 
for space-homogeneous sound propagation assuming a 
stochastic source of the form F{t)I{t) and including a 



2 



general causal form of Mori-Lee dissipation as well as 
linearizing in time-derivatives. The time integral of v(t) 
will then be the sound-level L{t) = log(7(t)//o) with a 
suitable reference level Iq expressed in decibel, which is 
an easily measurable quantity for traffic flows. With- 
out lack of generality the memory function T{t) and the 
stochastic background F(t) can be assumed to be sym- 
metric under time-reversal since v{—t) = —v{t). The 
fluctuations of F{t) are supposed to be stationary with a 
correlation compatible with the generalized fluctuation- 
dissipation theorem [T^ . i.e. 

POO 

{F{t + T)F{T)) =AfT{t) =Af dLupFiuj)cos{ujt) , 



(2) 

where ppii^) is the power spectrum of the noise 0. In 
terms of conventional Laplace transform techniques one 
now finds that 



M{f{t-t')Q{t-t') + f{t'-t)Q{t'-t)) 



Here 



I{t) 



1 

27ri 



7+ioo 



ds- 



r{s) 



(3) 



(4) 



where 7 is a sufficiently large positive number and where 
(t!^)o = (w(O)w(O)). The initial data leads in general to 
the normalization condition /(O) = 1 and stationarity of 
{v{t)v{t')) implies Af = {v'^)o as well as /(O) = 1. Eq.©, 
furthermore, leads to the following general equation for 
the fluctuations of the sound-level L{t): 



AL{tr^{{L{t)-L{0yf) = 



where we have defined the functions 



7+ioo . St 



ds e'^-l 



2m J^^ioc s .3 + r(s) 



(5) 



(6) 



/2W = TT- / -T— • (7) 



and 



For a stationary process the fluctuations in the sound- 
level L{t) is determined by the second term in Eq.lO 
with a unique short-time behavior, i.e. 



AL{t) = WJ\f+0{t') 



(8) 



and in principle the parameter J\f can therefore be deter- 
mined uniquely by a study of the sound fluctuations at 
sufRciently small time-scales. It is clear from Eqs.(|5ll and 
(0 that the large-time behavior of AL(t) is controlled by 
the properties of r(s) at small s -"-f r(s) ~ s"~^ at 
small s, i.e. T{t) ~ t~" at large t, one predicts that 



AL{t) ~ if a < 2. If e.g. r(s) = ^ for all s, with ^ 
real and positive, we obtain the conventional large-time 
behavior 



AL{t) 



(9) 



with a diffusion constant D = M/S,. If r(s) ~ C'^/s, 
with ^ real and positive, corresponding to a 1/uj power 
spectrum of pf{^) at small oj, we obtain 



AL{t) 



(10) 



apart from possible oscillatory and/or logarithmic cor- 
rections. If a < 2 we obtain from Eq.Q the asymp- 
totic limit (w^)oo = liint^oo(w(i)^) = A/". In the lim- 
iting case a = 2 we write r(s) = cs for small s, where 
c 7^ —1 is a constant. We then find that {v'^)oo = {{v^)o~ 
AA)/(l + c)2+AAand AL{ty = i{v^)o + cJV)tyil + c)\ 
For a stationary correlation {v{t)v{t')) we then see that 
M = {v^)oo{= {v^)o) and AL(t)2 = AAtV(l + c). In 
passing we notice that for a > 2, Eq.© leads to the 
asymptotic value (w^)oo = {v'^)o and according to Eq.© 
the asymptotic fluctuations are also independent of M in 
this case, i.e. AL(t)^ = {v'^)ot'^ for large values of t. In 
Ref. 13 it was also argued that in general A/" 7^ (i'^)oo, 
if a > 2, which was interpreted as an actual breakdown 
of the fluctuation-dissipation theorem. The results pre- 
sented above summarizes what we need in order to ana- 
lyze the large-time behavior of traffic flows. 



III. WAVELET ANALYSIS OF TRAFFIC FLOW 
DATA 



The new experimental method we suggest is to make 
use of the noise of sound as produced by traffic flows. The 
large-time behavior of the sound data obtained from such 
traffic flows can be characterized in terms of the scaling 
properties of suitably averaged wavelet amplitudes of an 
input signal S{t), i.e. 

W[S]{Xa) = X"W[S]{a) , (11) 

where H is the Hurst exponent and W [S] stands for the 
wavelet transform of the signal S{t) given in form of a 
time-series. Our normalization is such that H = 0.5 cor- 
responds to conventional Brownian motion. The averag- 
ing procedure we make use of follows the suggestion of 
Ref.[ll|, i.e. at a given scale a we average over the abso- 
lute values of the signal S{t) in the wavelet domain and 
hence M^[S'](a) ~ AS{t) with t ~ a. AL{t)^ then scales 
like t^^ at large times. If we, in a similar manner, de- 
fine a Hurst exponent Hp of the stochastic background 
F{t), i.e. {F{t + t)F(t)) scales like t^^^ at large times, 
then H = —Hp = a/2 according to Eqs.Q and |(SJ|. In 
our numerical work we have made use of the Daubechies 
wavelet family of order four |0, . 

As an example we have used data from a one-lane sim- 
ulation with identical cars located at sites, with an ini- 
tial occupancy probability p, of a lattice with 10* lattice 
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FIG. 1: This figure shows the typical scaling of averaged 
wavelet coeflicients for the \og{v{t) + 2) data as obtained from 
a one-lane simulation with a traffic intensity of 0.53 cars per 
unit time at the point of sound measurements with = 2^"^ 
data points. All identical cars moves in this case with the 
maximum speed of 5 lattice units per unit time. The density 
of cars per lattice site is p — 0.1016. For large time-scales we 
obtain an expected Hurst coefficient H — 0.5 corresponding 
to normal Einstein diffusion. 



points with periodic boundary conditions and N — 2^^ 
data points in time in order to illustrate our method. 
The actual value of cars per site, i.e. the car density per 
lattice site p, change very little during the simulation and 
approaches p in the limit of a large lattice. The cars have 
a hard core repulsion so that only one car at a time oc- 
cupies a given lattice site and a speed limit restricts the 
car velocities. In order to simulate the approach to an 
optimal traffic flow, a given car either speeds up or slows 
down depending on the distance to the car ahead. The 
velocity increases with one lattice unit per unit time if 
this is allowed by the distance to the car ahead. Other- 
wise the velocity is decreased with one lattice unit per 
unit time (for details of the implementation of such a 
simulation see e.g. Ref.j^). In order to make a concrete 
model we assume that the sound produced from a car 
at the site of the detection device leads to a sound-level 
L{t) = log(2 + v{t)), where v{t) is the velocity of the car 
at the site in lattice units per unit time. If there is no 
car present we put v{t) = —1, which simply corresponds 
to no sound at all. 

In Fig^we show the typical scaling of averaged wavelet 
coefficients for S{t) = \og{v{t) + 2) at a traffic intensity of 
0.53 cars per unit time. The maximum speed of the cars 
is 5 lattice units per unit time and p — 0.1. The actual 
number of cars per lattice site, p, is found to be 0.1016. 
For large time-scales we obtain an expected Hurst expo- 
nent H — a/2 — 0.5 corresponding to normal Einstein 
diffusion. It is interesting to observe that this scaling 
can be modelled by a constant power spectrum pf{u>), if 
< < LiJmax, with a cutoff parameter ujmax that still 



FIG. 2: This figure shows the scaling of wavelet coefficients as 
in Fig^but with a high traffic density p = 0.8064 per lattice 
site. The average velocity is in this case < v{t) >= 0.24 
lattice units per unit time. The traffic intensity is 0.1936 
cars per unit time. The wavelet analysis reveals two types of 
characteristic behavior. There is a range of a scale with an 
expected H — 0.5 corresponding to normal Einstein diffusion 
and a large-time scaling with H = 0.9 corresponding to a 
superdiffusive behavior. 



opens up a window for various finite-time behavior |l2l |. 
We intend to discuss this feature of our model in more 
detail elsewhere. Even though the first data points cor- 
responds to a linear t scaling we do not attribute much 
statistical significance to this behavior. In the next sim- 
ulation we increase the density of cars, i.e. we choose 
p = 0.8, but keep the rest of the parameters the same. 
Due to the presence of traffic jamming we expect an in- 
crease of the fluctuations of noise at larger time scales 
and therefore a large diffusion constant, i.e. a superdif- 
fusive behavior. Due to the larger car density one may 
expect a lower mean velocity and we actually find that 
< V >— 0.24 lattice units per unit time corresponding to 
a lower traffic intensity of 0.194 cars per unit time as com- 
pared to the simulation of FigQ] As presented in Fig[21 
we find, as expected, a normal diffusive behavior at small 
scales. At larger scales we find a Hurst exponent close to 
H ~ 1.0, i.e. a clear superdiffusive behavior. A model 
for such a superdiffusive behavior has been discussed in 
Ref.|l2l| in terms of a constant power spectrum pf{^) in 
a finite frequency interval < LUmm l£ l£ ^max with 
parameters which again opens up a window for various 
finite-time behavior. 

We have also performed a field measurement of a real 
traffic flow during a period of a high traffic intensity. The 
traffic data is based on a sample of = 2^^ measure- 
ments, one per second, on a two-lane finite road segment 
with opposite traffic but with a common crossroad. The 
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model for such a behavior a power spectrum of the form 
ppioj) — B ex.p{—uj/Q,)/Lu, vahd for u) > oJmin > 0, where 
ujjnin plays the role of an infrared cutoff and i? is a nor- 
malization constant. f2 >> Wmm is a corresponding ul- 
traviolet cutoff. The memory function is then of the form 
r{t) ~ B login exp{-C)/u}min{l + {nt)^/'^), valid for 
times t < exp(— C)/wmi„. Here C ~ 0.57721566... is the 
Euler-Mascheroni constant. One then finds that T{s) ~ 
Bujmin/ s^- The importance of a 1 //-noise in traffic flows 
has been noticed before (see e.g. Refs. [l9ll20| ). Here we 
presented a method in which such a behavior is easily 
revealed in sound fluctuations of traffic flows. 

IV. FINAL REMARKS 
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FIG. 3: This figure shows the scaling of wavelet coeflicients 
for the sound-level L{t) using sound data from a real two- 
lane traffic flow with A'^ = 2^^ measurements each with a 
time-window of one second. The conventional diffusive be- 
havior if = 0.5 is clearly seen at small scales. At large scales 
the sound-level becomes peaked around its mean value with a 
white-noise power spectrum for AF{t)^ which, formally, cor- 
responds to _ff = — 0.5. At intermediate scales we uncover 
a l//-noise which indicates a self-regulatory behavior of the 
trafiic flow. 



required instrument is a low-cost standard decibel meter 
with sound data automatically stored 0|. We expect 
conventional diffusion at short time-scales. At sufficiently 
large scales we expect a constant sound-level with a fiat 
power spectrum of the sound fluctuations since the pres- 
ence of the common traffic obstacle should average out in 
the data at large time scales. A fiat power spectrum for 
the fluctuations of the sound-level leads to AL{t)'^ ~ 1/t 
at large times. For medium scales it is not entirely clear 
in what manner the driver of the vehicles adapts to each 
other. As seen in FigOwe flnd the expected small scale 
H = 0.5 behavior. At large scales we flnd a Hurst expo- 
nent H — —0.5, i.e. formally r{t) ~ t at large times. This 
behavior corresponds to the expected flat, i.e. white- 
noise, power spectrum of the fluctuations AL{t)'^. It, 
however, appears that our model above breaks down for 
Hp = —H > 0, i.e. a < 0, and we do not consider 
this behavior in more detail. At intermediate scales wc 
find a low-frequency power spectrum pf{'^) — i.e. 
the famous 1// noise (see e.g. Ref.Q and references 
therein) , and the sound noise fluctuations are constant as 
in Ea. (|10|l . In order to be more specific, we consider as a 



In summary, we have considered some aspects of the 
large-time structure of traffic flows. We have suggested 
a cost-effective method to collect traffic data in terms 
of the noise of sound the traffics flows actually produce. 
The scaling behavior observed can easily be modelled in 
terms of a general form of stochastic differential equation 
with memory effects included which opens up an avenue 
of detailed analysis. In our presentation of the method 
retardation effects, Doppler shifts, frequency dependence 
of the sound and other features, which can be inferred 
from sound observations, have all been neglected but can 
be used in more realistic situations to extract details of 
the monitored vehicles, flow rates and other characteris- 
tics of traffic flows. Our data from a realistic traffic situ- 
ation also suggest the presence of a 1// noise, a feature 
which has been suggested elsewhere to be of great im- 
portance in the self-regulatory behavior of intense traffic 

flows m. 

Even though we have applied our methods in order 
to model properties of traffic flows, it appears that the 
considerations outlined in the present paper have a much 
broader range of applicability. We intend to come back 
to such issues elsewhere. 
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